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摘 要 : 湖 冰 物候 变化 特征 是 区 域 气候 变化 的 敏感 指示 器 之 一 。 近 几 十 年 来 ,由 于 全 球 变 暖 和 人 为 


活动 的 影响 ,中 亚 地 区 的 气候 发 生 了 显著 变化 ,准确 监测 湖 冰 物候 对 于 了 
的 科学 意义 。 通 过 对 中 亚 地 区 7 个 大 型 湖泊 (卡拉 库 尔 湖 、 巴 尔 喀 什 湖 、 成 海 、 阿 拉 湖 、 帝 又 


有 重 


解 中 亚 地 区 气候 变化 其 


油 、 查 蒂 尔 - 科 尔 湖 以 及 马 卡 科 尔 湖 ,面积 >100 km?) 2000 一 2020 年 的 长 期 地 表 反 射 率 数据 、 气 象 
数据 以 及 湖泊 资料 的 综合 分 析 , 利 用 GIS 相 关 技 术 探 讨 其 湖 冰 物 候 特 征 及 其 影响 因素 。 结 果 表 
明 :(1) 中 亚 地 区 的 湖泊 在 9 月 中 名 至 11 月 上 名 期 间 开 始 结 冰 ,11 月 底 到 12 月 底 会 完全 封冻 ,湖泊 
平均 冻结 时 间 为 35 d; 湖 冰 在 3 月 底 至 $ 月 中 开始 消融 ,4 月 上 至 6 月 上 会 完全 消融 ,湖泊 平均 消融 
时 间 为 18 d, (2) 2000 一 2020 年 中 亚 7 个 湖泊 中 有 5 个 湖泊 开始 冻结 日 期 呈现 延 后 的 趋势 ,平均 延 
后 速率 为 4.86 d- (10a) ,巴尔 喀什 湖 开 始 冻结 日 期 呈现 提前 趋势 ,提前 率 为 1.44 d.(10aj"。 完 全 消 
融 日 期 呈现 提前 的 趋势 ,平均 提前 率 为 2.90 d.(10a)'。7 个 湖泊 的 平均 湖 冰 存 在 期 为 171 d, 其 中 有 4 
个 湖泊 湖 冰 存在 期 呈 缩 短 趋 势 ,完全 冻结 期 呈现 整体 缩短 趋势 ,其 中 巴尔 喀什 湖 缩短 最 明显 ,缩短 速 


率 为 9.02 d.(10 a), (3) 中 亚 7 个 湖泊 湖 冰 的 冻结 -消融 空间 模式 主 
逐渐 冻结 ,消融 时 从 湖岸 到 对 岸 ;湖水 自 湖岸 冻结 至 


分 为 两 类 :湖水 从 两 岸 向 湖 心 
对 尾 , 越 先 冻结 的 湖区 越 先 消融 ,(4) 中 亚 地 区 湖 


泊 湖 冰 物 候 变 化 受到 湖泊 本 身 ( 海 拔 和 面积 ) 和 气候 (气温 和 降水 量 ) 等 多 个 因素 的 影响 ,气温 是 影响 


湖 冰 物候 的 关键 因素 ,气温 越 高 , 湖 冰 存在 期 越 短 ;面积 主 


影响 湖泊 的 冻结 日 期 ,面积 越 大 ,湖泊 的 


湖 冰 存在 期 越 短 ; 湖 冰 物 候 也 表现 出 一 定 的 海拔 依赖 性 , 随 着 海拔 的 升 高 ,湖泊 的 湖 冰 存 在 期 越 长 。 
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气候 变化 是 世界 范围 内 一 个 重要 的 共同 课 
题 ，“。 据 IPCC 第 六 次 报告 称 ,相对 于 1850 到 1900 
年 ,全球 地 表 温 度 已 经 升 高 1.09 °C ,并 呈现 出 加 速 
上 升 趋势 2 。 和 气候 变化 对 全 球 陆地 和 海洋 生态 系统 
产生 了 巨大 而 深远 的 影响 ,同时 给 人 们 生产 生活 带 
来 了 诸多 不 利 影 响 , 因 此 明确 气候 变化 响应 机 制 迫 
在 眉 睫 ”。 自 然 界 中 有 很 多 介质 可 以 用 来 指示 气候 
变化 ,比如 海平 面 . 冰 川 \ 冻 土 湖泊 等 ,湖泊 在 全 球 
水 循环 中 发 挥 着 至 关 重 要 的 作用 ,其 湖泊 面积 水 
位 以 及 湖 冰 物候 等 都 可 以 作为 区 域 气候 变化 的 指 
标 ', 其 中 湖 冰 物候 对 区 域 气 候 变 化 响应 异常 敏 
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PERO 。 高 纬度 地 区 的 湖泊 会 随 着 天 气 变 冷 而 产生 
湖 冰 ”", 湖 冰 物 候 变 化 对 气候 、 环 境 和 社会 经 济 活 
动 会 产生 重要 的 影响 。 此 外 , 湖 冰 物候 与 湖水 
能 量 平 衡 变化 有 直接 的 关系 “, 可 以 很 好 地 反映 区 
域 乃 至 全 球 气候 变化 ”““。 过 去 几 十 年 来 ,北半球 
的 湖 冰 正在 持续 减少 *, 而 且 随 着 未 来 气温 升 高 ,这 
种 情况 将 会 更 加 严重 "。 因 此 ,研究 湖 冰 物候 特征 
及 相关 因素 ,对 于 进一步 认识 湖泊 对 全 球 气 候 变 化 
的 响应 规律 和 反馈 机 理 ,具有 十 分 重要 的 科学 意义 。 

湖 冰 物候 研究 早期 主要 采用 野外 实地 观测 的 
方法 来 获取 湖 冰 物候 数据 。 但 这 种 方法 容易 受到 
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天 气 .环境 .技术 的 影响 ,无 法 获取 准确 .长 时 间 序 
列 大 范围 的 湖 冰 物候 数据 。 近 年 来 , 随 着 遥感 技 
术 的 鞍 勃 发 展 , 利 用 高 时 空 分 辩 率 的 遥感 影像 来 获 
取 湖 冰 物 候 信 息 已 经 成 为 一 种 研究 趋势 ”"。 基 于 
遥感 的 湖 冰 物候 监测 主要 有 微波 遥感 和 光学 遥 
感 。 微 波 遥 感 又 有 被 动 微 波 和 主动 微波 之 分 ,其 中 
被 动 微波 遥感 具有 高 时 间 分 辩 率 (每 天 2 次 或 更 
多 ) ,但 空间 分 辨 率 相 对 较 低 (>10 km) ,因此 被 动 微 
波 适 合 监测 大 型 湖泊 的 湖 冰 变化 ,例如 青藏 高 原 上 
的 青海 湖 “” ,以 及 加 拿 大 的 大 能 湖 和 大 奴 湖 ”和 
中 国 、 俄 罗斯 之 间 的 界 湖 兴 凯 湖 ”。 而 主动 微波 遥 
感 具 有 较 高 的 空间 分 辨 率 ,ERS( 方 位 方向 <30 m, 
距离 方向 <26.3 m) 被 用 于 监测 湖泊 结 冰 过 程 和 冰 
厚度 ~; Radarsat(1~100 m) 被 用 于 监测 湖泊 的 冻 
结 和 消融 过 程 汪 。 相 比 被 动 微 波 遥 感 , Eh 
感 的 时 间 分 辨 率 较 低 (ERS 为 3 d,Radarsat 为 24 d), 
不 能 满足 对 湖 冰 进行 高 频次 检测 的 需求 了”。 光 学 
遥感 如 MODIS 影 像 ,其 同时 具有 高 时 间 和 高 空间 分 
SES (1 d, 空 间 分 辨 率 为 250 m) ,被 广泛 应 用 于 监测 
湖 冰 物候 变化 。 如 :分 雪 楠 等 ?使 用 MODIS BACHE , 
分 析 了 2000 一 2020 年 色 林 错 湖 冰 物候 特征 及 其 影 
响 因 素 。Yao 等 外 利用 MODIS FI Landsat 数据 提取 
T 2000—2011 年 可 可 西里 地 区 22 个 湖泊 的 湖 冰 物 
候 , 并 对 相关 的 影响 因子 进行 了 讨论 。 吴 艳 红 等 ?” 
使 用 MODIS 数 据 ,提取 了 纳木错 湖 2000 一 2015 年 间 
的 湖 冰 物候 数据 ,基于 相关 模型 重建 了 1963 一 2018 
年 纳木错 湖 的 湖 冰 物 候 序 列 , 分 析 了 60 多 年 来 纳 木 
错 湖 冰 物 候 的 变化 特点 。 

先前 湖 冰 物 候 的 研究 多 集 中 在 高 海拔 上 纬度 
地 区 ,对 于 中 亚 地 区 湖 冰 物候 的 研究 较 少 。 中 亚 位 
于 欧 亚 大 陆 中 部 ,75% 的 地 区 属于 干旱 半 干 旱 的 大 
陆 性 气候 ,其 中 包括 数 千 个 湖泊 ,这些 湖泊 对 中 
亚 地 区 生态 系统 的 可 持续 性 和 人 类 福 社 有 着 重要 
的 意义 。 中 亚 湖 泊 的 相关 人 研究 主要 集中 在 湖泊 
面积 以 及 水 位 与 水 量 的 变化 。 如 Che 等 ”利用 ter- 
raPulse™ 月 度 Landsat 衍 生 的 地 表 水 域 范围 数据 集 和 
HydroLAKES 数据 集 , 对 2000 一 2015 年 中 亚 湖 泊 面 
积 的 时 空 变化 进行 了 研究 。Huang 等 “利用 SRTM 
和 Landsat 影 像 构 建 了 面积 -水 位 -水 量变 化 的 经 验 
模型 ,并 估计 9619 个 湖泊 中 最 大 水 量 和 最 小 水 量 。 
Hu 等 5 利用 3 个 GRACE 卫星 数据 集 和 5 个 全 球 水 
文 模型 对 2003 一 2014 年 中 亚 干旱 区 陆地 储 水 量变 


化 进行 研究 。 对 于 中 亚 地 区 湖泊 湖 冰 物候 的 研究 
相对 较 少 ,而 且 精 度 较 差 。 因 此 ,本 文 基于 MODIS 
地 表 反 射 率 数据 ,提取 了 2000 一 2020 年 中 亚 地 区 淹 
泊 的 湖 冰 物候 信息 ,分 析 了 湖 冰 物候 变化 特征 ,并 
利用 气象 数据 和 现 有 湖泊 资料 ,讨论 了 了 气温、 降水、 
湖泊 面积 以 及 海拔 对 湖 冰 物候 的 影响 。 


1 研究 区 概况 


中 亚 一 般 意 义 上 包括 土库曼 斯 坦 、 哈 萨克斯 
JH .吉尔吉斯 斯 坦 .塔吉克 斯 坦 和 乌 效 别克 斯 坦 。 
中 亚 的 地 势 总 体 上 东南 高 .西北 低 ,气候 干燥 ,降水 
稀少 ,昼夜 温差 较 大 ”*” , 占 世 界 干 旱 半 干旱 地 区 总 
面积 的 34%5 ,是 全 球 最 大 的 非 地 带 性 干旱 区 馈 
(图 1)。 中 亚 干 旱地 区 缺 水 严重 ,可 供 开 发 的 淡水 
资源 很 少 ” ,因此 对 于 中 亚 地 区 ,湖泊 是 非常 重要 
的 水 资源 后 。 中 亚 湖 泊 总 面积 超过 88000 km? ,其 中 
约 3000 个 湖泊 面积 >1 km’, 40 多 个 湖泊 面积 超过 
100 km’ ”。 由 于 独特 的 气候 ,中 亚 大 多 数 湖泊 的 
水 源 来 自 冰 川 融化 .山区 降水 和 河流 径流 汪 。 

其 中 选择 7 个 湖 冰 存在 期 稳定 的 并 且 面 积 > 
100 km 湖泊 ( 云 量 较 少 ,像素 误差 分 类 较 少 ,湖泊 并 
结 消融 明显 ) 作 为 研究 湖泊 ,它们 分 别 是 卡拉 库 尔 
湖 、 巴 尔 喀什 湖 、 咸 海 、 阿 拉 湖 、 帝 桑 泊 、 查 带 尔 - 科 
尔 湖 以 及 马 卡 科 尔 湖 。 表 1 为 7 个 湖泊 的 详细 信 
息 ,包括 海拔 高 度 面积 大 小 以 及 湖泊 类 型 。 


2 数据 与 方法 


2.1 数据 

2.1.1 MODIS 地 表 反 射 率 数据 MODIS 是 Terra 和 
Aqua 两 条 卫星 上 重要 的 传感器 ,本 文 所 使 用 是 MO- 
DIS 逐日 地 表 反 射 率 产 品 MODO09GQ (Terra) 和 
MODO9GQ( Aqua) ,该 数据 时 间 分 辨 率 为 1 d, 空 间 分 
辨 率 为 230 m。 利 用 RS 和 6GIS 技 术 对 遥感 数据 进行 
处 理 ,提取 7 个 湖泊 的 湖 冰 物候 。 

2.1.2 Landsat 数 据 Landsat 是 美国 航空 航天 局 与 
美国 地 质 调查 局 的 联合 计划 ,该 计划 提供 了 现 有 的 
世界 上 最 长 的 连续 空 基 观测 记录 ,从 计划 开始 到 现 
在 , 共 发 射 了 9 颗 卫 星 。 本 研究 用 到 的 Landsat 数据 
(时 间 分 辨 率 为 16 d, 空 间 分 辨 率 为 30 m) 主 要 包 
括 :Landsat7、Landsat8( 表 2) ,通过 获取 MODIS 和 
Landsat 同时 过 境 的 数据 ,来 检验 MODIS 数 据 提取 湖 
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图 1 研究 区 湖泊 分 布 
Fig. 1 Distribution of lakes in the study area 


R1 研究 中 使 用 的 湖泊 信息 


Tab.1 Lake information used in the study 


湖泊 名 称 地 理 位 置 面积 /km 年 份 海拔 /m 湖泊 类 型 
卡拉 库 尔 湖 39°02’24"N ,73?25' 12"E 380 2015 3914 咸水湖 
巴尔 喀什 湖 46?10'27"N ,74°20'25"E 16996 2021 342 东 为 咸 西 为 淡 
查 带 尔 - 科 尔 湖 40937'25"N ,75°18'20"E 181 2020 3530 咸水湖 
马 卡 科 尔 湖 48^45'23"N ,85°45'29"E 455 2019 1445 咸水湖 
阿拉 湖 46?10'17"N,81?35'15"E 2650 2016 347 咸水湖 
毅 桑 泊 48°00'07"N , 84?00'05"E 1810 2021 388 淡水 湖 
咸 海 45953'25"N ,60?23'21"E 3300 2008 39 咸水湖 

表 2 研究 中 使 用 的 Landsat 数 据 
Tab.2 Landsat data used in the study 

陆地 卫星 影像 日 期 (年 -月 -日 ) 湖泊 名 称 陆地 卫星 影像 日 期 (年 -月 -日 ) 湖泊 名 称 

Landsat7 2016-11-17 卡拉 库 尔 湖 Landsat8 2018-12-09 马 卡 科 尔 湖 

Landsat7 2016-12-03 巴尔 喀什 湖 Landsat8 2018-12-25 查 带 尔 - 科 尔 湖 

Landsat7 2017-01-10 查 带 尔 - 科 尔 湖 Landsat8 2019-05-02 威海 

Landsat7 2017-05-12 马 卡 科 尔 湖 Landsat8 2019-05-18 毅 桑 泊 

Landsat7 2017-05-28 毅 桑 泊 Landsat8 2019-11-26 巴尔 喀什 湖 

Landsat7 2017-11-20 咸 海 Landsat8 2019-12-12 阿拉 湖 

Landsat7 2017-12-06 阿拉 湖 Landsat8 2019-12-28 卡拉 库 尔 湖 

Landsat7 2017-12-22 卡拉 库 尔 湖 Landsat8 2020-01-29 毅 桑 泊 

Landsat8 2018-01-07 巴尔 喀什 湖 Landsat8 2020-05-04 成 海 

Landsat8 2018-05-15 阿拉 湖 Landsat8 2020-05-20 卡拉 库 尔 湖 

冰 物 侯 属 性 的 精度 。 平均 气温 与 降水 量 , 在 中 亚 气 候 研 究 中 表现 出 良好 
2.1.3 和 再 分 析 气 象 数据 由 于 中 亚 地 区 气象 站 稀少 。 的 适用 性 ”*”。 本 研究 获取 了 2000 一 2020 年 的 气 


且 连 续 观 测 数据 周期 短 ,气象 站 数据 不 足以 文 撑 研 
究 湖 冰 物 候 变 化 的 因果 分 析 。 因 此 ,本 研究 使 用 再 
分 析 气 象 数 据 (CRUTS v4.05) 来 分 析 湖 冰 物 候 变 化 
的 影响 因素 。CRUTS v4.05 数 据 ( 时 间 分 辨 率 为 1 个 
月 ,空间 分 辨 率 为 0.$") ,包括 1901 一 2020 年 逐 月 的 


温 以 及 降水 数据 来 对 湖 冰 物候 变化 进行 相关 性 
AT TS 

22 方法 

2.2.1 湖 冰 信息 提取 冰 在 可 见 光 与 近 红 外 波段 的 
反射 率 高 ,而 水 在 可 见 光 和 近 红 外 波段 的 反射 率 较 


564 f 44k UR 47% 


低 ” ,因此 根据 冰 和 水 在 红 光 与 近 红 外 波段 的 反射 
率 差 异 ,利用 直方 图 和 目 视 解 译 ,找到 合适 闪 值 来 
区 分 冰 和 水 ,该 方法 为 装 值 法 ”"。 具 体 公 式 为 : 


. . {Band1 - Band2 >a 
Lake ice = i e (1) 


式 中 :Bandl、Band2 分 别 为 MODI 影像 250 m 的 红 
光 、 近 红外 波段 ;ab 为 浆 值 , 当 Bandl-Band2>a、 
Band1>b 时 , 则 可 以 认为 是 湖 冰 。 

NDSI ( Normalized difference snow index) 77 1X: th, 
可 以 用 于 提取 湖 冰 信 息 ”。 本 文 使 用 NDSI 方 法 提 
取 Landsat 数 据 的 湖 冰 特征 信息 ,其 公式 如 下 : 


. Band2 - Band5 
DIDSES Band2 + Band5 (2) 


式 中 :Band2 为 绿 光 波段 ;Band5 为 短波 红外 波段 。 
依据 计算 出 来 的 NDSI 结 合 反 射 率 直方 图 以 及 目 视 
解 译 最 终 确 定 浆 值 , 进 而 区 分 出 湖 冰 和 湖水 。 
2.2.2 湖 冰 时 间 属 性 定义 湖 冰 时 间 属 性 具体 分 为 : 
湖泊 开始 冻结 日 期 (FUS) 完全 冻结 日 期 (FUE) 、 开 
台 消融 日 期 (BUS) 完全 消融 日 期 (BUE)。 湖 冰冻 
结 期 (FUD ) 指 湖泊 从 开始 冻结 直至 完全 冻结 所 用 
的 时 间 ; 湖 冰 消 融 期 (BUD) 指 湖泊 从 开始 消融 直至 
完全 消融 所 用 的 时 间 ;湖泊 湖 冰 存 在 期 (ICD) 指 从 
湖泊 开始 冻结 到 完全 消融 这 上 段 时 间 ; 湖 泊 完 全 冻结 
期 (CFD) 指 从 湖泊 完全 冻结 至 开始 消融 的 这 段 时 
期 。 一 般 情况 下 ,在 湖泊 冻结 期 间 ,湖泊 边缘 附 
近 会 出 现 部 分 未 冻结 的 区 域 , 在 湖 冰 消融 过 程 中 ， 
湖岸 处 也 可 能 会 存在 零星 的 湖 冰 “3” 。 此 外 ,由 于 
天 气 变 化 湖泊 边界 不 匹配 、 像 素 分 类 错误 以 及 不 
可 避免 的 噪声 (尤其 是 云层 ) 导 致 的 反复 冻结 和 
消融 ,会 降低 湖 冰 信息 提取 的 准确 性 。 因 此 ,本 研 
究 根 据 Kropáek 等 “3 的 方法 ,将 湖 冰 占 湖泊 面积 的 
10% 和 90% 作 为 闷 值 进行 湖 冰 物候 的 提取 。 具 体 公 
式 如 下 : 


湖 冰 物 候 提 取 = 
FUS， 湖 冰 面 积 二 湖泊 面积 的 10% 
FUE， 湖 冰 面 积 = 湖 泊 面 积 的 90% (3) 
BUS， 湖 冰 面 积 =< 湖 泊 面 积 的 90% 
BUF， 湖 冰 面 积 和 湖泊 面积 的 10% 


3 结果 与 分 析 


3.1 MODIS 与 Landsat 提 取 湖 冰 物 候 信 息 的 交叉 
验证 
为 了 评估 阅 值 法 提取 MODIS 湖 冰 信 息 的 精度 ， 


本 文 使 用 了 Landsat 数 据 进行 了 验证 。 在 本 人 研究 中 ， 
选用 了 测量 平均 绝对 误差 (MAE) VJ E FE TE CR?) Fl 
偏差 (Bias) 等 指标 来 评估 了 MODIS 数 据 提取 湖 冰 物 
候 的 精度 。 首 先 , 选 取 了 20 幅 没有 云层 干扰 以 及 湖 
泊 处 于 冻结 与 消融 过 程 中 的 Landsat 影 像 ,然后 分 别 
计算 出 每 一 幅 影 像 中 湖 冰 面积 占 整 个 湖泊 面积 的 
百分比 。 最 后 将 这 些 比 例 与 MODIS 数据 提取 的 湖 
冰 履 盖 百 分 比 进行 对 比 。 

图 2 为 Landsat 数 据 与 MODIS 数据 提取 湖 冰 面 
积 比 例 的 对 比 , 其 结果 为 : 尼 为 0.999,MAE 为 0.92%， 
Bias X 1.1596 , 这 表示 MODIS 数据 根据 阔 值 法 提取 
中 亚 湖泊 的 湖 冰 面积 精度 较 高 ,可 用 于 提取 中 亚 地 
区 大 型 湖泊 的 湖 冰 物候 信息 。 图 3 为 基于 NDSI 与 
闵 值 法 提取 2017 年 12 月 22 日 卡拉 库 尔 湖 的 湖 冰 
言 息 ,其 中 Landsat 影 像 使 用 NDSI 方 法 提取 的 湖 冰 
面积 为 297.69 km’, MODIS 数据 基于 阅 值 法 提取 的 
湖 冰 面积 为 286.15 km? ,2 种 提取 方法 的 误差 约 为 
0.58%。 
3.2 2000 一 2020 年 中 亚 7 湖 的 湖 冰 物 候 特 征 分 析 

基于 MODIS 影像 ,使 用 闪 值 法 结合 人 工 目 视 解 
译 , 得 到 2000 一 2020 年 中 亚 地 区 所 选 湖泊 的 湖 冰 物 
候 平 均 特 征 ( 表 3、 图 4)。 

从 图 4 和 表 3 可 以 看 出 ,7 个 湖泊 从 每 年 9 一 11 
月 开始 结 冰 , 其 中 阿拉 湖 和 查 带 尔 - 科 尔 湖 开始 冻 
结 日 期 是 早 于 7 个 湖泊 的 平均 冻结 日 期 。 阿 拉 湖 开 
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& 6 Bias=1.15% 
B n-20 
AS 
AE 
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MODIS 湖 冰 面 积 /% 
IE: MAE 为 平均 绝对 误差 ;R 为 测定 校正 ; 


Bias 为 债 差 刀 为 影像 幅 数 。 下 同 。 
图 2 Landsat 数 据 与 MODIS 数据 湖 冰 面积 比例 对 比 
Fig. 2 Comparison of lake ice area ratio between Landsat 
data and MODIS data 
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(a) Landsat (2017-12-22) (b) MODIS (2017-12-22) 
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图 3 NDSI-S BUE Epok A B. 
Fig. 3 Extraction of lake ice information by NDSI and threshold method 


表 3 2000 一 2020 年 中 亚 地 区 所 选 湖泊 平均 湖 冰 物 候 统计 


Tab.3 Average lake ice phenology of selected lakes in Central Asia from 2000 to 2020 /d 
湖泊 名 称 开始 冻结 ”完全 冻结 ”开始 消融 ”完全 消融 冻结 期 消融 期 ” 湖 冰 存在 期 完全 冻结 期 

巴尔 喀什 湖 333 360 453 471 29 19 140 94 
查 带 尔 - 科 尔 湖 283 330 499 518 49 20 237 170 
马 卡 科 尔 湖 324 336 498 505 13 8 182 163 
阿拉 湖 261 334 450 482 75 33 223 117 
卡拉 库 尔 湖 315 359 493 516 45 23 202 136 
咸 海 336 - 365 463 - - 127 - 

毅 桑 泊 326 335 471 479 10 8 153 137 
平均 值 311 342 477 491 35 18 171 126 


注 :采用 年 积 日 法 ,第 333 d 为 当年 11 月 27 日 ,第 453 d 为 次 年 3 月 28 日 。 


Or @ p 00 始 冻结 日 期 最 早 ,威海 开始 冻结 日 期 最 晚 。7 个 湖 
te Au Gd URN IS. FUS PUE — 泊 的 平均 完全 冻结 时 间 是 在 当年 的 11 月 底 到 12 月 
500 - 平均 BUE 


底 ,平均 冻结 期 为 35 d. 
7 个 湖泊 的 湖 冰 在 次 年 的 3 一 4 月 开始 消融 ,其 

中 巴尔 喀什 湖 和 阿拉 湖 日 期 是 早 于 7 个 湖泊 的 平均 
so} 消融 日 期 。 阿 拉 湖 日 期 开始 消融 时 间 最 早 , 而 查 带 
ee 尔 - 科 尔 湖 开始 消融 的 时 间 最 晚 。7 个 湖泊 的 平均 完 

0 er 全 消融 时 间 是 在 次 年 的 46 月 ,平均 消融 期 为 18 


ee a S SS d。7 个 湖泊 中 ,有 6 个 湖泊 的 冻结 期 比 消融 期 长 ,说 
2 明 中 亚 地 区 所 选 湖泊 冻结 的 速度 比 消融 的 速度 慢 。 
T 湖泊 名 称 


7 个 湖泊 在 湖 冰 存在 期 和 完全 冻结 期 上 差异 很 
注 :FUS 为 开始 冻结 日 期 ;FUE 为 完全 冻结 日 期 ;BUS 为 开始 消融 日 大 。 湖 泊 平 均 湖 冰 存在 期 为 171 d, 其 中 咸 海 湖 冰 的 
HH. Mp, Im qi 期 。 采 积 日 法 ,第 Mx, MZ m — ds os . NS 
c pres ct 1127 湖 冰 存在 期 最 短 为 127 d, 查 蒂 尔 - 科 尔 湖 冰 存 在 其 
E (人 o H o 
E :湖泊 平均 完 冻结 其 ; 其 
图 4 2000 一 2020 年 中 亚 地 区 所 选 湖泊 平均 湖 冰 物候 特征 TEPER 26 Bese 126 LATA 
We AK BERI 56 2 2 TR 170 d, 巴 尔 喀什 淹 


Fig. 4 Average lake ice phenological characteritics of 
selected lakes in Central Asia from 2000 to 2020 完全 冻结 期 最 短 为 94 do 
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3.3 2000 一 2020 年 中 亚 7 湖 的 湖 冰 物候 变化 趋势 

2000 一 2020 年 7 个 湖泊 湖 冰 物候 变化 趋势 存在 
较 大 的 差异 ( 表 4、 图 5~6)。 从 表 4 和 图 5 中 可 以 看 
出 ,2000 一 2020 年 除了 巴尔 喀什 湖 [1.44 d+ (10a) ] 
和 卡拉 库 尔 湖 [0.30 d - (103)! | AY FUS = [8 X ni t 
势 外 , 剩 下 5 个 湖 的 冻结 日 期 都 是 呈 延 后 趋势 ,平均 
延 后 率 为 4.86 d (10a) ', FCP Ar aie 2K BL EI Bix 
为 明显 ,推迟 率 为 18.00 d- (103), BUE BR T KALE 
尔 湖 外 缓慢 延 后 外 ,其 余 的 呈现 提前 的 趋势 ,平均 
提前 率 为 2.90 d(10a)', 其 中 阿拉 湖 提前 最 明显 , 提 
前 率 为 5.70 d- (10a), 

2000 一 2020 年 查 带 尔 - 科 尔 湖 、 马 卡 科 尔 湖 、 阿 
拉 湖 和 斋 桑 油 的 湖 冰 存在 期 呈现 出 缩短 趋势 。 例 
如 ,阿拉 湖 湖 冰 存 在 期 缩短 趋势 为 8.7 de (103) ; 5 
卡 科 尔 湖 、 阿 拉 湖 以 及 背 桑 泊 的 完全 冻结 期 和 湖 
冰 存 在 期 都 呈现 缩短 趋势 。 在 7 个 湖泊 中 ,除了 
阿拉 湖 和 威海 的 冻结 期 与 消融 期 的 趋势 相反 ,其 
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余 5 个 湖泊 在 冻结 期 和 消融 期 方面 表现 出 相同 的 
趋势 。 


4 讨论 


4.1 湖泊 冻结 -消融 空间 模式 特征 
通过 对 中 亚 地 区 所 选 6 个 湖泊 的 冻结 -消融 过 
程 统 计 分 析 , 湖泊 冻结 -消融 空间 格局 分 为 2 种 : 
(1) 冻结 自 湖岸 一 侧 延 伸 至 另 一 侧 , 越 先 冻结 的 湖 
区 越 先 消融 。 符 合 这 种 空间 格局 的 湖泊 包括 :卡拉 
库 尔 湖 .巴尔 喀什 湖 、 阿 拉 湖 . 毅 桑 泊 以 及 马 卡 科 尔 
湖 。 为 了 直观 展示 这 种 模式 的 湖 冰 物候 演变 过 程 ， 
本 人 研究 选用 巴尔 喀什 湖 为 例 , 如 图 7 和 图 8 所 示 ， 
由 于 巴尔 喀什 湖西 部 淡水 区 的 盐 度 通常 低 于 0.5 g: 
L ,而 东部 半 咸 水 区 的 盐 度 可 以 高 达 3 g LOME, 
所 以 巴尔 喀什 湖 于 11 月 下 旬 左 右 先 在 湖西 南方 癌 
开始 冻结 , 湖 冰 从 西南 向 东北 缓慢 覆盖 ,到 12 月 下 


表 4 2000—2020 中 亚 地 区 所 选 湖泊 湖 冰 物候 变化 趋势 


Tab.4 Change trend of lake ice phenology in selected lakes in Central Asia from 2000 to 2020 


湖泊 名 称 开始 冻结 完全 冻结 开始 消融 

巴尔 喀什 湖 -1.44 7.01” -2.01 
查 蒂 尔 - 科 尔 湖 18.00” -0.40 2.40 
马 卡 科 尔 湖 0.70 0.40 -4.90 
阿拉 湖 3.00 6.00" -0.90 
卡拉 库 尔 湖 -0.30 6.80™ -1.50 
威海 1.70 - - 

毅 桑 泊 1.00 0.90 -4.40 


注 :*、**、*** 表 示 通 过 P 值 为 0.1、0.05 ,0.01 的 显著 性 检验 。 


gx 
v 威海 
图 例 
| 湖 冰 物候 变 化 率 /d*(10a)? 
本 ? 


mm 7459 


/d+(10a)" 


完全 消融 冻结 消融 期 ”完全 冻结 期 ” 湖 冰 存在 期 
-0.58 8.33" 1.44 -9.02 0.87 
-1.20 -19.00™ -5.00 4.10™ -19.20 
-5.20' -0.40 -0.30 -5.30' -5.90" 
-5.10" 3.00 -4.70" -7.00" -8.70' 

0.70 7.007 2.20 -8.30 1.007 
-1.50 - - - -3.20 
-4.00* -0.10 -0.40 -4.60' -5.00 


”巴尔 喀什 湖 


注 : 湖 冰 物候 变化 率 柱状 图 朝 上 为 延 后 趋势 , 朝 下 为 提前 趋势 。 
图 5 2000 一 2020 中 亚 地 区 所 选 湖泊 湖 冰 物 候 变 化 趋势 空间 分 布 
Fig. 5 Spatial distribution of lake ice phenology change trend in selected lakes in Central Asia from 2000 to 2020 
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图 6 2000 一 2020 年 中 亚 所 选 湖泊 湖 冰 物候 年 际 变化 


Fig.6 Interannual variation of lake ice phenology in selected lakes in Central Asia from 2000 to 2020 


A FCB A i ,整个 湖 冰 冻结 期 比较 漫长 。 湖 亲 
消融 时 ,主要 从 湖 的 西南 岸 先 消融 ,然后 向 东北 逐 
渐 消 融 , 到 4 月 中 旬 , 湖 冰 完 全 消融 (2) 湖水 从 两 
岸 向 湖 心 逐渐 冻结 ,消融 时 从 一 侧 到 另 一 侧 ,这 种 
模式 的 湖泊 有 查 蒂 尔 - 科 尔 湖 。 图 9 为 查 蒂 尔 - 科 
尔 湖 湖 冰 物候 演化 过 程 , 查 带 尔 - 科 尔 湖 于 10 月 上 
旬 左 右 先 在 湖南 北 两 岸 形成 岸 冰 , 湖 六 由 两 岸 向 湖 
心 逐 渐 昌 延 ,到 12 月 下 旬 完 全 被 冰 上 履 善 ,整个 湖 冰 
冻结 期 比较 长 。 湖 冰 消 融 时 ,主要 从 湖 的 南岸 先 消 
融 ,由 南 向 北 迅 速 消融 ,到 5 月 中 名 , 湖 冰 完全 


消融 。 
42 气象 因素 对 湖 冰 物 候 变 化 的 影响 

湖 冰 物候 的 形成 受到 气候 条 件 和 湖泊 理化 特 
征 的 影响 “5 ,其 中 对 湖 冰 物候 产生 影响 的 气象 因 
素 有 :气温 降水 太阳 辐射 . 风速, 风 轩 等, 以往 在 
其 他 地 区 的 研究 发 现 , 气 温 是 影响 湖 冰 物候 的 关键 
因素 中。 为 了 研究 中 亚 地 区 选 定 湖泊 的 湖 冰 物候 
变化 与 气候 变化 之 间 的 响应 规律 ,本 研究 分 析 了 气 
温和 降水 量 与 湖 冰 物候 之 间 的 相关 性 。 

为 了 验证 再 分 析 气 象 数 据 的 准确 性 ,本 研究 使 
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图 7 2010—2011 年 巴尔 喀什 湖 冻 结 -消融 过 程 
Fig.7 Freezing-thawing process of Balkhash Lake from 2010 to 2011 
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图 例 
巴尔 喀什 湖 盐 度 /gL 

0 100 km 
i >; L 


图 8 巴尔 喀什 湖 盐 度 分 布 
Fig.8 Salinity distribution map of Balkhash Lake 


用 2013 年 6 月 到 2020 年 12 月 卡拉 库 尔 湖 气象 站 的 “中 亚 地 区 的 气候 变化 状况 。 

月 平均 气温 与 CRUTS v4.05 对 应 的 数据 进行 了 交叉 图 11 为 湖 冰 物候 与 气温 、 降 水 的 相关 性 。 由 图 
验证 ,如 图 10 所 示 尼 为 0961,MAE 为 0211% Bias 。 11 可 知 , 湖 冰 的 存在 期 完全 冻结 期 .开始 和 完全 融 
为 1.807% ,这 表明 CRUTS v4.05 数 据 集 可 用 于 讨论 ”化 期 与 年 平均 气温 有 较 强 的 相关 性 。 年 平均 气温 
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图 9 2010—2011 年 查 蒂 尔 - 科 尔 湖 冻结 -消融 过 程 
Fig. 9 Freezing-thawing process of Chatir Kol from 2010 to 2011 
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CRU 和 气温 /"C 
注 :n 为 对 应 数据 的 个 数 。 
图 10 CRUTS v4.05 数 据 集 精度 验证 


y=0.976x+2.138 


Fig. 10 Precision validation of CRUTS v4.05 dataset 
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较 低 时 , 湖 冰 的 存在 期 和 完全 冻结 期 较 长 ,年 平均 
气温 较 高 时 , 湖 冰 的 存在 期 较 短 。 此 外 ,中 亚 湖泊 
开始 和 完全 融化 的 日 期 与 年 平均 气温 呈 较 强 的 负 
相关 。 这 表明 气温 是 中 亚 湖泊 湖 冰 物 候 变化 的 决 
定 因素 。 此 外 ,分 析 结 果 还 表明 ,降水 对 中 亚 地 区 
湖泊 湖 冰 物候 的 变化 影响 不 大 。 

4.3 湖水 理化 性 质 对 湖 冰 物候 变化 的 影响 

43.1 湖 冰 物 候 将 征 与 湖泊 面积 之 间 的 关系 湖泊 
面积 会 影响 湖 冰 物候 的 变化 ”。 网 12 为 近 20 a7 个 
湖泊 的 平均 面积 与 FUS、BUE ICD 的 相关 性 。 从 图 
中 可 以 看 出 ,中 亚 地 区 7 个 湖泊 面积 与 FUS 呈 正 相 
关 , 与 BUE 呈 负 相 关 , 这 表明 湖泊 面积 越 大 ,湖水 开 
始 结 冰 的 时 间 愈 迟 ,融化 愈 快 , 湖 冰 的 存在 期 就 


降水 量 


0.0 NONEM 0.1 00 -00 


FUS FUE BUS BUE FUD BUD CFD ICD FUS FUE BUS BUE FUD BUD CFD ICD 
湖 冰 物候 属性 


注 :FUD 为 冻结 期 ;BUD 为 消融 期 ;CFD 为 完全 冻结 期 ;ICD 为 湖 冰 存 在 期 ;负数 代表 负 相 关 。 
图 11 湖 冰 变 化 物候 与 气象 因子 变化 的 相关 性 热 图 


Fig. 11 Correlation heat map of lake ice change phenology and changes in meteorological factors 
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图 12 中 亚 地 区 湖泊 湖 冰 物候 与 湖泊 面积 关系 图 


Fig. 12 Relationship between lake ice phenology and lake area in Central Asia 


图 13 为 2000 一 2020 年 中 亚 湖泊 面积 与 FUS 变 
化 率 、BUE 变化 率 、ICD 变化 率 的 相关 性 ,从 图 中 可 
以 看 出 面积 与 FUS 变 化 率 、BUE 变化 率 、ICD 变化 率 
的 相关 系数 分 别 为 -0.42 .0.24、-0.45 , 湖泊 面积 后 
FUS 和 1CD 变化 率 呈 现 负 相关 ,这 说 明湖 油 面 积 越 
大 ,湖泊 的 FUS 和 ICD 的 变化 速率 越 慢 ;湖泊 面积 三 
BUE 变化 率 呈 现 正 相关 ,表明 湖 油 面积 越 大 ,BUE 
变化 速率 越 快 。 

43.2 湖 冰 物候 特征 与 湖泊 海拔 之 间 的 关系 湖面 


5j FUS, BUE ICD 的 相关 性 。 从 网 中 可 以 看 出 海拔 
与 FUS 星 负 相 关 , 与 BUE 和 ICD 呈 显著 正 相 关 。 这 
说 明湖 面 海拔 越 高 ,湖水 开始 结 冰 的 时 间 越 早 , 消 
融 的 时 间 越 晚 , 湖 冰 存在 期 越 长 。 

图 15 为 2000 一 2020 年 中 亚 湖泊 海拔 与 FUS 变 
化 率 、BUE ZERA ICD 变化 率 的 相关 性 ,从 图 中 可 
以 看 出 海拔 与 FUS、BUE、ICD 的 相关 系数 分 别 为 
0.49 .0.63 .0.43 ,都 呈现 出 正 相 关 趋 势 , 这 说 明 随 着 
海拔 的 升 高 ,湖泊 的 FUS、BUE 以 及 ICD 的 变化 速率 


海拔 会 影响 湖 冰 物候 的 变化 ”。 图 14 为 湖泊 海拔 ”都 在 加 快 。 
21 | (a) 面积 与 FUS 变 化 率 .(b) 面积 与 BUE 变 化 率 6r (© 面积 与 ICD 变 化 率 


1 
0 


al E aa 
m 15 E " = 2 
a = . = 
S 12 x ~) p 入 0 
Mt 9 ^ a> AM =} 
wd = -3 y-9.364x—3.089 à . y=-3.153x-3.586 
& 6r.. . y=—4.407x+5.497 E 4 : R2-0.059 A m QU. R2-0.206 
D 3 0.. R0174 巴 O -6 E 
= Ofer . ELS d m 一 SF u cor 
n =i 
0 5000 10000 15000 20000 0 5000 10000 15000 20000 0 5000 10000 15000 20000 
面积 /ko 面积 (km 面积 /ko 
图 13 中 亚 地 区 湖泊 湖 冰 物 修平 均 变 化 率 与 湖 当面 积 关 系 图 
Fig. 13 Relationship between the average change rate of lake ice phenology and lake area in Central Asia 
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图 14 中 亚 地 区 湖泊 湖 冰 物 候 与 湖泊 海拔 关系 图 
Fig.14 Relationship between lake ice phenology and lake altitude in Central Asia 
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图 15 中 亚 地 区 湖泊 湖 冰 物 修平 均 变 化 率 与 湖泊 海拔 关系 图 
Fig. 15 Relationship between the average change rate of lake ice phenology and lake altitude in Central Asia 


4.4.4 无 量 纲 化 分 析 

无 量 纲 化 是 一 种 数据 分 析 方 法 , 旨 在 将 数据 集 
中 的 多 个 变量 转换 为 一 组 无 单位 的 变量 ,以 便 更 好 
地 理解 和 分 析 数 据 。 这 个 过 程 可 以 减少 数据 的 复 
AAS VE ,提高 数据 的 处 理 效 率 , 并 提高 数据 分 析 的 精 
度 和 准确 性 。 本 文采 用 标准 化 进行 无 量 纲 分 析 , 具 


(4) 


WP: x H n 个 训练 样本 中 第 j 个 特征 值 组 成 的 问 
量 ; jy 为 训练 样本 中 的 均值 ; o 为 训练 样本 的 标 
准 差 。 

441 单个 湖泊 的 无 量 岗 化 分 析 以 巴尔 喀什 湖 为 
例 , 将 湖 冰 物 候 信息 、 冬 半年 气温 、 降 水 量 、 面 积 采 
用 标准 化 方式 进行 无 量 纲 化 处 理 , 然 后 进行 多 元 回 
归 分 析 , 进 而 分 析出 冬 半 年 气温 、 降 水 量 、 面 积 对 湖 
冰 物 候 的 影响 程度 。 具 体 分 析 结 果 如 下 : 

FUS -0.457X' -0.030Y + 0.247Z" + 8.734E - 11 

FUE =0.657X"" ~ 0.037Y + 0.427Z" -7.511E - 11 
BUS = -0.109X + 0.150Y - 0.013Z + 6.560E - 11 
BUE- -0.306X + 0.013Y - 0.0682 - 2072E-12 (5) 
FUD -0.285X^ - 0.023Y + 0.237Z' +3.584E - 11 

BUD = -0.170X^ - 0.237Y + 0.099Z - 2.867E - 11 


CFD = -0.503X" + 0.022Y - 0.160Z —3.986E - 11 
ICD 2 -0.517X^ * 0.114Y - 0.302Z + 2.648E - 11 


式 中 :XX 为 冬 半年 气温 ;了 为 降水 量 ;Z 为 面积 ;* 、**、 
**x* 表 示 通 过 P 值 为 0.1、0.05、0.01 的 显著 性 检验 。 
结果 表明 ,气温 是 影响 湖 冰 物候 的 关键 因素 ,面积 
主要 影响 湖 冰 物 候 的 FUE 和 CFD, 降 水 量 对 湖 冰 物 
候 的 影响 相对 较 小 。 

4.4.2 不 同 湖 泊 的 无 量 岗 化 分 析 以 6 个 湖泊 
2000 一 2020 年 的 平均 湖 冰 物候 .平均 湖泊 面积 以 及 
湖泊 海拔 进行 无 量 纲 处 理 , 然后 进行 多 元 回归 分 
Vr ,进而 分 析出 面积 和 海拔 对 湖 冰 物候 的 影响 程 
度 。 具 体 分 析 结 果 如 下 : 


FUS 2 0.357X - 0.199Y + 3.324F - 11 

FUE = 0.204X - 0.578Y - 9.637E - 11 

BUS = -0.660X + 0.224Y - 3.737E - 11 
BUE--1.158X" - 0.914Y' - 1.524E - 10 (6) 
FUD - -0.220X + 0.346Y - 5.774E - 11 

BUD =0.376X -0.617Y + 6.389E - 11 

CFD = -0.254X + 0.306Y + 2.177E — 10 

ICD 7 -0.047X + 0.520Y - 8.661E - 11 


式 中 :了 为 海拔 ;了 为 面积 。 结 果 表 明 , 除 了 FUE 主要 
受 湖泊 面积 影响 外 ,其 他 湖 冰 物候 主要 受 海 拔 影响 。 
4.5 现 有 研究 结果 对 比 

通过 文献 搜索 查阅 发 现 , 人 研究 中 亚 地 区 湖 冰 物 
候 变 化 的 文 草 非 常 少 ,我 们 统计 了 现 有 研究 结果 和 
本 研究 的 结果 的 差异 ( 表 5)。Hou 等 “使 用 了 MO- 
DIS 地 表 温 度数 据 ,通过 设置 一 定 国 值 来 提取 湖 冰 


X5 本 研究 与 现 有 研究 结果 对 比 


Tab. 5 Comparison of the results of this study with existing research 


作者 人 研究 区 域 人 研究 数据 数据 分 辨 率 BEYE El% 
Hou" ”北半球 中 MODIS 地 表 温 。 1d;lkm 20 
纬度 地 区 上 度 产品 80 
本 文 中 亚 五 国 “MODIS 地 表 反 。 1d;250m 10 
射 率 产 品 90 


时 间 跨 度 优点 缺点 
2002 一 2020 年 ”一 天 可 以 使 用 四 幅 影 像 ” "DG ETE SIS BO [HC 
进行 湖 冰 属 性 的 提取 设 定 ,提取 精度 较 差 
2000 一 2020 年 ”数据 空间 分 辩 更 高 , 湖 “采用 动态 阔 值 提取 
冰 属 性 提取 精度 更 高 湖 冰 信息 , 耗 时 长 


5/2 


的 范围 ;而 本 研究 采用 MODIS 地 表 反 射 率 数据 , 通 
过 设置 动态 国 值 来 提取 湖 冰 。 前 者 虽然 可 以 一 天 
获取 多 景 数据 来 减少 云 的 影响 ,但 由 于 不 同 湖泊 目 
号 理化 性 质 (如 盐 度 等 ) 的 差异 ,使 用 固定 国 值 融 提 
取 湖 冰 物 候 属 性 会 存在 一 定 的 误差 。 而 本 研究 采 
用 250m 的 反射 率 数 据 , 具 有 更 高 的 空间 分 辨 率 , 且 
SUAS YRC. BML, HAS BEBE Tes o 


5 结论 


AS CHEF MODIS 5 Landsat 数据 并 结合 相关 的 
气象 数据 ,分析 了 中 亚 地 区 7 个 湖泊 的 长 期 湖 冰 物 
候 特征 及 其 影响 因素 , 主要 结论 如 下 : 

(1) 2000 一 2020 年 中 亚 的 湖泊 从 每 年 9 一 11 月 
开始 结 冰 , 湖 冰 在 3 一 6 月 逐渐 完全 消融 。7 个 湖 整 
体 平 均 冻 结 期 和 平均 消融 期 为 33 4 和 18 d, WYKE 
均 完 全 冻结 期 和 平均 存在 期 为 126d4 和 171 do 

(2) 2000 一 2020 年 中 亚 7 个 湖泊 中 有 5 个 湖泊 
开始 冻结 日 期 呈现 延 后 的 趋势 ,平均 延 后 率 为 
4.86 d- (10a)'。 总 体 上 ,完全 消融 日 期 有 提前 的 趋 
A ,平均 提前 率 为 2.90 d*(10a)!'!。 其 中 有 4 个 湖泊 湖 
冰 存 在 期 呈 缩 短 趋势 , 平 均 缩短 率 为 9.90 d - (10a) '. 
完全 冻结 期 呈 整 体 缩短 趋势 ,平均 缩短 率 为 5.00 d- 
(10a)". 

(3) 中 亚 7 个 湖泊 湖 冰 的 冻结 -消融 空间 模式 
主要 分 为 两 类 :湖水 自 湖岸 冻结 至 对 岸 , 越 先 冻 结 
的 湖区 越 先 消融 ;湖水 从 两 岸 向 湖 心 逐渐 冻结 , 消 
RAEN Gi Fe EDO Fe o 

(4) 中 亚 地 区 湖 冰 物候 变化 受到 多 种 因素 的 综 
合影 响 ,包括 湖泊 自身 的 特征 (如 海拔 .面积 等 ) 和 
气候 变化 (如 气温 和 降水 量 等 )。 气 温 是 影响 湖 冰 
物候 的 关键 因素 ,气温 越 高 , 湖 冰 存在 期 越 短 ; 面积 
越 大 ,湖泊 的 湖 冰 存在 期 越 短 ; 随 着 海拔 的 升 高 , 湖 
泊 的 湖 冰 存在 期 越 长 。 
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Temporal and spatial variations of lake ice phenology in large 
lakes of Central Asia from 2000 to 2020 


ZHAO Mingjie", WANG Ninglian", SHI Chenlie", HOU Jingqi” 

(1. Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Northwest University, Xi an 710127, 
Shaanxi, China; 2. Institute of Earth Surface System and Hazards, College of Urban and Environmental Sciences, Northwest 
University, Xi an 710127, Shaanxi, China; 3. State Key Laboratory of Tibetan Plateau Earth System, Environment and 
Resources (TPESER), Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing. 100101, China) 


Abstract: The phenology of lake ice is a sensitive indicator of regional climate change. Through comprehensive 
analysis of long-term surface reflectance data, meteorological data, and lake information for seven large lakes 
(Karakul Lake, Balkhash Lake, Aral Sea, Alakol Lake, Zaysan Lake, Chatir Kol Lake, and Markakol Lake) with 
an area greater than 100 km’ in Central Asia from 2000 to 2020, GIS-related technologies were used to explore 
the characteristics of lake ice phenology and its influencing factors. The results are as follows: (1) Lakes in Cen- 
tral Asia began to freeze from mid-September to early November and completely froze from late November to 
late December, with an average freezing time of 35 days; lake ice began to melt from late March to mid-May and 
would completely melt from early April to early June, with an average melting time of 18 days. (2) From 2000 to 
2020, the start dates of ice formation in five of the seven lakes in Central Asia exhibited a delayed trend, with an 
average delay rate of 4.86 days per decade, whereas the start date of Balkhash Lake exhibited an advancing trend, 
with an advancing rate of 1.44 days per decade. The analysis suggests that this may be due to a decrease in the an- 
nual average temperature in the winter half of the year. The complete melting dates showed an advancing trend, 
with an average advancement rate of 2.90 days per decade. The average ice-covered period for the seven lakes 
was 171 days, with four of the lakes exhibiting a trend of shortening of the ice- covered period. The complete 
freezing period shows an overall trend of shortening, with Balkhash Lake exhibiting the most significant reduc- 
tion, with a rate of 9.02 days per decade. (3) The spatial pattern of the formation and melting of lake ice in the 
seven lakes in Central Asia can be mainly divided into two categories: the lake water gradually freezes from both 
sides to the center and melts from the lake shore to the opposite side, or the lake water freezes from the shore to 
the opposite side and the earlier freezing, lake area melts the sooner. (4) The lake ice phenology changes in Cen- 
tral Asia are influenced by multiple factors such as lake characteristics (altitude and area) and climate (tempera- 
ture and precipitation). Temperature is the key factor affecting lake ice phenology, and the higher the temperature, 
the shorter the ice-covered period. The area primarily affects the freezing date of the lake, and the larger the area, 
the shorter the ice-covered period. As the altitude increases, the ice-covered period of the lake extends. 


Key words: lake ice phenology; climate change; Central Asia; MODIS 


